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Abstract 

We present ab-initio density functional (DFT) calculations of the vi- 
brational spectra of neutral Magnesium phthalocyanine (MgPc) molecule 
and of its Raman scattering intensities. 

1 Introduction 

Organic materials are attractive for electronic industry, because they promise 
large versatility based on organic chemistry, low price based on very common 
elements used as building blocks and interesting mechanical properties. There is 
quite a vigorous research activity in organic materials for electrical applications 
- from insulators trough semiconductors to conductors and even superconduc- 
tors. Here we are concerned with metal phthalocyanines (MPc). Electronic 
structure calculations ^ and recent experimental observation of the possibility 
to dope these materials by adding electrons |21 resembles somewhat the situ- 
ation of fuUerenes. In our previous paper we speculated about possible phase 
diagrams for electron doped MgPc |3] . While these possibilities are under active 
experimental considerations, we noticed that for the commonest experimental 
diagnostic, namely Raman scattering there is no reference calculation of either 
modes nor Raman intensities (with exception of ZnPc |14|). In this paper we 
focused our interest on the vibrational properties of MgPc, which we wish to 
study ab-initio. While our final aim will be to compute spectra of doped and 
undoped molecules, this work will be restricted to neutral undoped MgPc, which 
has not been studied so far. 



2 Technical details. 

All our calculations were done using the PWscf software package 0], which is 
a plane-wave basis set DFT implementation. We used the LDA approximation 
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with Slater approximation for exchange and Perdew-Zunger functional for cor- 
relations effects. In order to lower the plane wave energy cut-off we used the 
non-local RRK J3 ultra-soft pseudopotentials [S| . The kinetic energy cut-off was 
35Ry for the wave functions basis set and 280Ry for the charge density basis set. 
The plane-wave basis set assumed periodic boundary conditions. The unit cell 
had dimensions 21x21xllA, therefore including enough vacuum to represent 
the isolated molecule. 

By energy optimization we found the equilibrium molecular structure with 
structural parameters as in Table [3 The molecule has symmetry group 
and is planar. The electronic structure agrees well with former calculations of 
Liao et. al.pj. Some levels near the HOMO-LUMO gap are listed in the table 
□ 

We calculated the vibrational spectra of neutral, undoped MgPc by means of 
the density functional perturbation theory . The dynamical matrix was calcu- 
lated in Cartesian coordinates. In the isolated molecule with A'' atoms there are 
3A'^ — 6 genuine vibrations while 6 modes corresponds to translation-rotational 
degrees of freedom with zero frequency. In reality due to periodic images of 
molecules there is a weak virtual interaction which can cause "libration" of the 
molecule at nonzero, even imaginary frequency. We eliminated these modes 
by transforming the dynamical matrix to internal coordinates and setting the 
corresponding dynamical matrix elements to zero. We checked that after this 
procedure we obtained eigenvectors of correct symmetry. Assignment of irre- 
ducible representation to all eigenvectors was done by projection on symmetry 
adapted bases of all linearly independent atomic displacement. For N — 57 
atoms there are 165 vibrations 

r„b = UAig + 13A2g + UBig + UB2g + 13Eg+6Aiu+8A2u + 7Biu + 7B2u+28Eu 

where Eg and i?„ are two-fold degenerated modes. The calculated frequencies 
are summarized in Fig. |3 

The Raman intensity is related to the change of molecular polarizability due 
to the deformation introduced by the vibration. We calculated the polarizability 
derivatives with respect to deformations in the static limit. 

With each vibrational mode qi is associated a Raman tensor given by 



r = (1) 



Here ^, v are Cartesian indices, a^^ is the polarizability tensor defined as 
= , U is the total energy and E an external electric field. Therefore 

one has to calculate third derivatives of total energy ^dE dxr "^^^ energy 
derivatives with respect to atomic displacements are forces acting on atomic 
nuclei, and are calculated via the Hellman-Feynman theorem. We extracted the 
force dependence on electric field by numerical differences. Applied fields were 
±1 X 10~'^, ±2 X 10"'^, ±4 X 10"'^ a.u.. Details about apphcation of electric 
field within periodic boundary conditions are the same as in Ref. '7'. There 
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are alternative methods for the Raman tensor calculation using linear response 
theory [SUHj- However numerical derivative method with finite field is computa- 
tionaly less demanding for the present case, where moreover the implementation 
is very simple. 

From the Raman tensor Eq.([3l we calculated for each mode the gas phase 
(angle averaged) Raman scattering cross section following |1(J[ [TT] in the form: 

dcTj 

In 

«^ 

Uq is frequency of incident light, is equilibrium occupation number for the 
initial vibrational state i at the given temperature T. The a' and 7'^ are 
isotropic and anisotropic parts of the Raman tensor. Both are invariant under 
rotations. Formula (j^J is valid for the most common experimental setup when 
incident ray, observed ray, and incident light polarization of electric field are 
perpendicular to each other. 

The Raman tensor of a particular mode must belong to the same irreducible 
representation as that vibration mode. In a coordinate system connected to 
the molecule as indicated on Fig. Q all symmetric tensors of second rank are 
divided according to irreducible representations in the form: 

/ a \ /cOO\ /OdO\ /OOe 

^ig : a , Big : -c , Bag : d , Bg : / 
\0 b J \000/ \000/ \e/0 

(3) 

Traceless irreducible representations possess only the anisotropic part, while the 
Aig mode has both parts. Weights in front of a' and 7' in formula |3 (45 and 
7 respectively) can be changed by observation geometry, thus providing a way 
to (partially) assign symmetry to vibration eigenmodes (see for example Ref. 

We note that for the Eg representation it is possible to define a single 
independent parameter yc^+Z^ because of the degeneracy and orthogonality 
of degenerated modes. 

3 Results and discussion 

In Table 121 we give our LDA optimized structure of MgPc. In comparison with 
data in references we see a reasonable agreement. The difference can be 

attributed to different density functional (GGA in Ref. ^). LDA is known to 
slightly overbind in covalent bonds. Kohn-Sham orbitals energies (Tab. are 
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'2 , T^,/^ 



(27ri^o)^ h{n, + 1) 45a' ■ + 7j 



45 



(2) 



(T* +T'' +T' 

V XX ' yy z. 



{{TL - Tlyf + {Tl, - Tl^Y + {T;^ - TJ,)' + 6[(T^,)' + (T^,)' + (t;,)']} 



exp 



kT 



given relative to vacuum zero taken as the SCF electrostatic potential far from 
the molecule. They are also in good agreement with Ref. ^ as far as we can 
read off values from fig. 2 of Ref. . 

Vibrational spectra of various metal-phthalocyanines were experimentally 
measured ^3] by means of Raman spectroscopy. In the case of ZnPc there is 
also a DFT calculation of the vibrational and Raman spectra In general 
the most intense peak was found to be in the interval 1500-1550cm~^ depending 
on the central atom. 

Our calculated vibration spectrum of MgPc is on Fig.||2J). Above the fre- 
quency line for Raman active modes we write the corresponding non-zero com- 
ponent of Raman tensor as defined by Eq.©. The simulated Raman spectrum 
according to formula |2] is shown on Fig. (PJ . To simulate spectrum for each 
intensity and frequency we added Gaussian with the spread a=5cm.~^ . We con- 
sidered the low temperature limit {rn = 0) and low incident light frequency 
limit. Our most intense Raman scattering mode is at 1587cm~l with symmetry 
Big. This is similar to ZnPc which is at 1516cm~^ with symmetry Bi. 

The results of the presented calculation should be compared with Raman 
data of MgPc, as soon as available. It will also serve as a starting point for 
future work on the electron doped molecules 
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level KS-eigenvalue [eV] symmetry 

LUMO+1 ^2^241 

LUMO -3.692 eg 

HOMO -5.132 aiu 

HOMO-1 -6.092 b2g 

HOMO-2 -6.446 62™ 

HOMO-3,-4 -6.526 e„ 

Table 1: Kohn-Sham eigenvalues spectra and irreducible representation of wave- 
functions near the HOMO-LUMO gap. Note that Cg and e„ orbitals are two-fold 
degenerated. 



bond DFT Ref. [J Experiment Kef. Hil this work 



Mg - Ni 


2.008 


1.990 


1.991 




1.377 


1.386 


1.359 


N2-C1 


1.335 




1.317 


Ci ~ C2 


1.465 


1.411 


1.445 


C2 — C2 


1.415 


1.468 


1.400 


C2 — 


1.395 


1.400 


1.379 


C3 — C4 


1.397 


1.399 


1.381 


C4 — C4 


1.406 


1.412 


1.392 


C-H 


1.090 


1.121 


1.092 




109.7 


109.5 


109.9 




127.5 


125.9 


127.1 


ON1-C1-C2 


108.6 


108.9 


108.4 



Table 2: Structure of MgPc. Bond lengths are in A, angles in degrees. Notations 
corresponds to labels on Fig. ^ 
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Figure 1: Schematic structure of the MgPc molecule. Orientation of coordinate 
system is important for definition of the Raman tensors by Eq. Ipj|l. 
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Figure 3: Calculated Raman spectrum of MgPc for very low excitation frequency 
and low temperature. 
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